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ABSTRACT 

Context. Modem radio telescopes allow us to record a large number of spectral channels. The application of a Fourier transform to 
spectropolarimetric data in radio continuum, Faraday rotation measure (RM) synthesis, yields the "Faraday spectrum", which hosts 
valuable information about the magneto-ionic medium along the line of sight. 

Aims. We investigate whether the method of wavelet-based RM synthesis can help us to identify structures of regular and turbulent 
magnetic fields in extended magnetized objects, such as galaxies and galaxy clusters. 

Methods. The analysis of spectropolarimetric radio observations of multi- scale targets calls for a corresponding mathematical tech- 
nique. Wavelets allow us to reformulate the RM synthesis method in a scale-dependent way and to visualize the data as a function of 
Faraday depth and scale. 

Results. We present observational tests to recognize magnetic field structures. A region with a regular magnetic field generates a 
broad "disk" in Faraday space, with two "horns" when the distribution of cosmic-ray electrons is broader than that of the thermal 
electrons. Each field reversal generates one asymmetric "horn" on top of the "disk". A region with a turbulent field can be recognized 
as a "Faraday forest" of many components. These tests are applied to the spectral ranges of various synthesis radio telescopes. We 
argue that the ratio of maximum to minimum wavelengths determines the range of scales that can be identified in Faraday space. 
Conclusions. A reliable recognition of magnetic field structures in spiral galaxies or galaxy clusters requires the analysis of data 
cubes in position-position-Faraday depth space ("PPF cubes"), observed over a wide and continuous frequency range, allowing the 
recognition of a wide range of scales as well as high resolution in Faraday space. The planned Square Kilometre Array (SKA) will 
fulfill this condition and will be close to representing a perfect "Faraday telescope". The combination of data from the Low Frequency 
Array (LOFAR, at low frequencies) and the Expanded Very Large Array (EVLA, at high frequencies) appears to be a promising ap- 
proach for the recognition of magnetic structures on all scales. The addition of data at intermediate frequencies from the Westerbork 
Synthesis Radio Telescope (WSRT) or the Giant Meterwave Radio Telescope (GMRT) would fill the gap between the LOFAR and 
EVLA frequency ranges. The Global Magneto-Ionic Medium Survey (GMIMS), planned with several single-dish telescopes at low 
angular resolution, will also provide good scale recognition and high resolution in Faraday space. 

Key words. Methods: observations - techniques: polarimetric - galaxies: magnetic fields - galaxies: spiral- galaxies: clusters: intra- 
cluster medium - radio continuum: galaxies 

1 ■ Introduction where B\\ is the line-of-sight magnetic field component measured 

in yuG, Hq is the thermal electron density measured in cm"^, and 
Modem radio telescopes are equipped with digital correla- the integral is taken along the line of sight through the region 
tors that allow us to record a large number of spectral chan- containing both magnetic fields and thermal electrons, where 



nels. While radio spectroscopy in total intensity is well- is measured in parsecs. Our adopted convention is that is pos- 

^ developed, the possibilities of spectropolarimetry in radio con- i^^y^ ^^en B is pointing towards the observer. 

tinuum have been exploredfor only a few years. The fundamen- Bum ( 1966) showed that the complex polarized intensity P 

tals were presented by |Bum| (ll966|), while the first application ^t different wavelengths A can be calculated from via a 

to multi-channel polarization data (data cubes) was presented by Fourier transform 

iBrentjens & de Br uyn (2005). ^co 

Faraday rotation measure (RM) synthesis generates the P(A^) = F(0)e^^'^'^ d0, (2) 

"Faraday dispersion function" or, in short, the "Faraday spec- ^'"^ 

trum" F((p), which is the (complex-valued) polarized intensity which means that the Faraday spectrum F is the Fourier trans- 
spectrum as a function of "Faraday depth" form of the complex polarized intensity: 



1 r 

0(x) = O.81 B\\(x)n^(x)dx\ (1) ^ 

^0 Fnn^itinn ^ 



P(A^)Q-^''^'^'dA^. (3) 



Equation © is the basis of RM synthesis. A major limitation 
emerges from the fact that P can be measured only for A^ > 
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Table 1. Spectral ranges of various radio telescopes and param- 
eters crucial for RM synthesis (see text for details). 



Telescope 


A 
m 


AA^ 


1^01 

rad/m^ 


|A0maxl 

rad/m^ 


(■^max/ 
■^min) 


LOFAR HBA 


1.25-2.73 


5.9 


0.59 


2.8 


4.8 


WSRT 


0.17-0.23 












+ 0.77-0.97 


0.91^ 


3.8 


110 


33 


GMRT 


0.21-0.30 
+ 0.47-0.52 












+ 0.87-0.98 


0.92^ 


3.8 


71 


22 


DRAO, Parkes, 


0.17-0.23 










Effelsberg 


+ 0.33-1.0 


0.97 


3.6 


110 


35 


(GMIMS) 












Parkes (S-PASS) 


0.12-0.14 


0.004 


870 


220 


1.4 


Arecibo 


0.20-0.24 


0.021 


165 


79 


1.4 


(GALFACTS) 












EVLA 


0.025-0.30 


0.089 


39 


5000 


144 


ATCA 


0.03-0.27 


0.072 


48 


3500 


81 


ASKAP 


0.21-0.27 


0.026 


130 


71 


1.6 


(POSSUM) 
(POSSUM 


0.21-0.42 


0.14 


25 


71 


4.0 


+ FLASH) 












SKA phase 1 


0.10-4.3 


18 


0.19 


310 


1850 


SKA phase 2 


0.03-4.3 


18 


0.19 


3500 


20500 



^ High sidelobes in Faraday spectrum expected owing to the large gaps 
in wavelength coverage 



The quantity F can be used to determine the Faraday depth 
and the intrinsic polarization angle of each component in the 
Faraday spectrum. As in classical spectroscopy, the interpreta- 
tion of this spectrum is not straightforward. In particular, there 
is no simple relation between Faraday depth and geometrical 
depth. Furthermore, the Faraday spectrum suffers from sidelobes 
of the main components caused by limited coverage of the wave- 
length space, and in the case of point-like "sources" these side- 
lobes can be removed by "RM clean" ( Heald.2QQ 9). 

RM synthesis is charac terized by three basic parameters 
(iBrentiens & de Bruvnll2005h : 

- the resolution S(p in Faraday space, which is inversely pro- 
portional to the coverage AA^ in wavelength (A^) space; 

- the maximum observable |0maxl of point-like sources in 
Faraday space, which is inversely proportional to the width 
of a single frequency channel; 

- the maximum width |A0niaxl of extended structures in 
Faraday space (Faraday-rotating and synchrotron-emitting 
regions), which is inversely proportional to the square of the 
minimum observation wavelength. Wide-band observations 
at long wavelengths yield high resolution in Faraday space 
but cannot detect extended structures. 

We argue below (Sect. 3) that the ratio of maximum to mini- 
mum wavelengths is another important parameter in helping us 
to recognize a range of scales in Faraday space. 

A region of any extent in physical space that is synchrotron- 
emitting but not Faraday-rotating (i.e. no thermal electrons) gen- 
erates a "point source" in Faraday space at = 0. A non-emitting 
region with thermal electrons and magnetic fields in front of the 
emitting region is called "Faraday screen"; it shifts the point 
source in Faraday space to a non-zero value of 0. 

RM synthesis has been applie d to spectropolarimetr i c data 
from the W esterbork telescope (de Bruvn & Brentiens 2005; 
Schnitzeler et al 2009; Hea ldetal. 2009; Maoetal. 2010; 
BrentiensI lioTll IPizzoet al.l [20111) . from the Parkes tele- 
scope and the Australia Telescope Compact Array (ATCA) 
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Fig. 1. Several typical examples of a region with and without re- 
versals of the regular magnetic field: the solid black line shows 
the Gaussian magnetic field distribution, the dot-dashed blue line 
stands for one reversal and the dashed red line shows the distri- 
bution with two reversals. Top: Magnetic field profile B(x) along 
the line of sight. Bottom: Faraday depth distribution 0(x) along 
the line of sight. 



Model 




hc/h^ 


Remark 


1 


V2 


1/V2 


Equipartition between ric and B 


2 


1 


1/2 


Similar distributions of ric and /ith 


3 


1/V2 


1/(2 V2) 


f7th more extended 


4 


2 


1 


He more extended 



Table 2. Models of the density distributions of cosmic-ray elec- 
trons (^c) and thermal electrons (^th) considered in this paper, he 
is the Gaussian scale-height (or scale-radius) of the distribution 
of cosmic-ray electrons, he that of the thermal electrons, and /zb 
that of the regular magnetic field B (see Section 4 for details). 



(iFeain et al.l l2QQ9l: iHarvev-Smith et al.l I201QI : IQ ' Sullivan etaP 
12012), from theVLA (I Van Eck et al.1 '201 1) and from the EVLA 
( Heesen et al. | [2QTTh . The AA^ coverages of these observations 
were small, so that the resolution in Faraday space was lim- 
ited. Neve rtheless, many Fa raday spectra have revealed complex 
structures (lBrentjensll201lh . By means of the application of RM 
synthesis, the fine- scale structure o f the magnetic field in the 
Milky Way around a local HI bubble (|Wolleben et al.ll2QlQah and 
around a supernova remnant (iHarvey-Smith et al. 12010^ could be 
measured. 

On the other hand, as RM synthesis is regularly ap- 
plied in radio astronomy, the number of questions discussed 
about its technical problems (algorithms, software) and prac- 
tical limitations for existing instruments similarly increases 
( Heald 2009"; 'Famsworth et al."2011^, 'Fri ck et al.ll201lt iLi et al.l 
2011; Andrecut et al. 2011; Bell et al. 20111). A combination of 
RM synthesis and two-dimensional image synthesis into three- 
dimensional (3-D) "Faraday synthesis" of data from synthesis 
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Fig. 2. Several typical examples of the Faraday spectrum in 
Faraday depth space F((p). Top: Gaussian magnetic field distri- 
bution for four diff'erent models (see Table 2): solid green line 
(model 1); dashed black line (model 2); dashed red line (model 
3); dot-dashed blue line (model 4). Bottom: The case of mag- 
netic field reversals and model 1): one reversal (dot-dashed blue 
line); two reversals (dashed red line). 



telescopes has been proposed as a possible way of improving 
the technique (Bell & EnBlin 2012). 

In this paper, we reconsider the problem of observations of 
large magnetized objects, such as galaxies or the intergalactic 
medium in galaxy clusters. Two factors are crucial in this con- 
text: the objects are very extended and composed of numer- 
ous structures of very diff'erent spatial scales, namely galaxies 
containing a disk, halo, spiral arms, turbulent star-forming re- 
gions, and supernova remnants, and galaxy clusters containing 
radio galaxies, jets, outflows, turbulent gas, and "relic" shock 
fronts. The multi-scale structure calls for the use of a mathe- 
matical technique developed in the framework of wavelet anal- 
ysis. Wavelets do not only provide a mo re robust algorithm on 
which to base the RM synthesis method (iFrick et al.ll201Ql ), but 
can also be used as an illustrative tool for the presentation of re- 
sults. We will determine here the information that can be recov- 
ered from RM synthesis based on wide-band observations with 
existing and future radio telescopes. 

For simplicity, we neglect in this paper the efl'ects of difl'erent 
angular resolutions and, in the case of synthesis telescopes, of 
difl'erent distributions of baselines (uv coverage) when combin- 
ing data from difl'erent telescopes. We also neglect the spectral 
variations in the polarized emission and the instrumental noise. 



2. Present and future radio telescopes 

Spectropolarimetric data can be obtained with the following 
present and future radio telescopes: 

- Low Frequency Array (LOFAR, Europe), which is the first 
online-connected synthesis telescope at low frequencies and op- 




-10 10 

(rad m 

Fig. 3. Model example of a turbulent magnetic field in the inter- 
galactic medium of a galaxy cluster. Panels from above: mag- 
netic field profile along the line of sight B(x), Faraday depth dis- 
tribution along the line of sight 0(x), and Faraday spectrum in 
Faraday depth space F((p). 



erates since 2011 fe tappers et al.ll20l Tl). Polarization from pul- 
sars has been detected in the "lowband" (10-80 MHz) and in the 
lower "highband" (1 10-190 MHz). Observations in the lowband, 
lower highband, and higher highband (about 170-240 MHz) 
need diff'erent setups and cannot be performed simultaneously. 
The lower highband has maximum sensitivity and is used to 
search for diffuse polarized emis sion, which is one task of the 
"magnetism key science project" (lBeckll20fQh . The lowband has 
lower sensitivity and signals suff'er from strong Faraday depolar- 
ization. In this paper, the full frequency range of the "highband" 
(HBA) of 110-240 MHz (or about 1.25-2.7 m) is considered. 

- Westerbork Synthesis Radio Telescope (WSRT, the 
Netherlands), which has successfully supplied spectropolarimet- 
ric data in the frequency range 1300-1763 MHz of a sample of 
nearby galaxies, t he SINGS survey ( Heald et al.l l2009), clusters 



(Pizzoeta 



( M ao et al 



1201 Ih . and to measure RMs of background sources 
20101) . The low resolution in Faraday space only 



allowed measurement of one component at most locations, 
except for the central regions of several galaxies where three 
components were detected. The lower frequency range 310- 
390 MHz (or 0.77-0.97 m) allows higher resolution in Faraday 
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space and hence is of great interest to RM syn thesis, e.g. as ap- 
plied to data of the diffuse Galactic foreground (Schnitzeler et al. 
2QQ9I) an d of galaxy clus ters (de Bruvn & Brentiens 2005; 
BrentiensI [201 ll: IPizzo et al.l 1201 ID . A survey of nearby galax- 
ies is underway in this frequency band. After installation of the 
APERTIF system, the WSRT will only observe around 1 GHz. 

- Giant Meterwave Radio Telescope (GMRT, India), a syn- 
thesis telescope that operates in the frequency ranges 150-156, 
236-244, 305-345, 580-640, and 1000- 1450 MHz. The polar- 
ization calibration of the GMRT is more difficult than that of the 
WSRT, and only a few successful detections of polarized signals 
from bright radio sources in the two hi ghest frequency bands 
have been reported so far (iJoshi & Chen galur 2010). In the fol- 
lowing we assume that the frequency range 305-1450 MHz can 
be used for future polarization observations, which is similar to 
that of the WSRT. 

- DRAG 26 m single dish (Penticton, Canada), for which a 
northern sky polarization surve y (GMIMS) is ongoing i n the fre- 
quency range 1277-1762 MHz (IWoUeben et al.ll2010bh . 

- Parkes 64 m single dish (Australia), where a southern 
sky survey (S-PASS) has bee n performed in the frequency 
range 2180-2420 MHz (ICarrettill2Q10l) . The southern part of the 
GMIMS survey at 300-900 MHz is planned. 

- Effelsberg 100 m single dish (Germany), at which a north- 
ern sky polarization survey (GMIMS) is planned in the fre- 
quency range 300-900 MHz. 

- Arecibo 305 m single dish (USA), where a deep spec- 
tropolarimetric survey (GALFACTS) of the sky visible from 
Arecibo is ongoing in the frequency range 1225-1525 MHz 
dGeorge et al. 2011). 

- Expanded VLA synthesis telescope (EVLA, USA), re- 
cently renamed as Jansky VLA, whose new correlator (WIDAR) 
allows measurements in the frequency range 1-12 GHz (0.025- 
0.30 m) with continuous coverage at four frequency settings (L, 
S, C, and X-band). Polarization data from the frequency bands 
at even higher frequencies are needed only for sources with 
|A0| > 5000rad/m^. The first polarization observations of the 
dwarf irregular galaxy ICIO in C band (4.5-7.8 G Hz) were still 
limited in resolution in Faraday space (Heesenet al. 2011). The 
full wavelength coverage available with the EVLA still needs to 
be investigated. 

- Australia Telescope Compact Array (ATCA), whose new 
correlator (CABB) allows measurements in the frequency range 
1.1-10 GHz (0.03-0.27 m) with an almost continuous cover- 
age at two frequency settings (L-hS and C-hX band). Complex 
Faraday spe ctra were observed in active galactic nuclei in the 
L-hS bands ( 0' Sullivan e'tani2012h . 

- Australian Square Kilometre Array Pathfinder (ASKAP), 
which is a synthesis telescope under construction in Western 
Australia (I Johnston et al.l [2QQ8I) . The all-sky polarization sur- 
vey POSSUM in the frequency range 1 130-1430MHz (or 0.21- 
0.265 m) is planned in 2013, together with surveys in HI and 
continuum. A combination of these data with other data from 
the all-sky transient survey FLASH in the frequency range 700- 
1000 MHz is under discussion; this would extend the frequency 
coverage to 700-1430MHz (or 0.21-0.428m) and improve the 
resolution in Faraday space by a factor of about five compared 
to data from the POSSUM survey alone. 

- Square Kilometre Array (SKA): Construction should start 
in 2016 in South Africa, with extensions into central Africa, and 
in Australia. The frequency coverage of the receiving systems 
under investigation for "phase 1" (which will be available from 
about 2018) is about 70 MHz-3 GHz (or 0.1-4.3 m) with contin- 
uous coverage, to be extended in "phase 2" to 10 GHz, opening 



excellent possibilities for RM synthesis. Measuring a dense grid 
of RM measurements is one of the key science projects for the 
SKA ( Gaensler et al. 2004; Beck 2010, 2011). 

Table [T] summarizes the properties of the present-day and fu- 
ture radio telescopes. The highest resolution in Faraday space 
(largest AA^) is achieved by LOFAR and the SKA, while the 
ATCA, ELVA, and the SKA provide the largest range of scales 
in Faraday space (largest (/Imax/^min)^, see Sect. [3]). 

The following combinations of data from synthesis tele- 
scopes yield an excellent frequency coverage and high angular 
resolution, and are considered in this paper: 

- LOFAR -h EVLA data of high angular resolution over a 
huge frequency range (110-240 MHz -h 1-12 GHz or 0.025- 
0.30 m -h 1.25-2.7 m), but the frequency gap hampers detection 
of certain structures in Faraday space. 

- LOFAR -h WSRT/GMRT -h EVLA data help to partly fill 
the gap in the data from LOFAR -h EVLA. Faraday structures 
can be detected on a large range of scales, similar to the planned 
capability of the SKA. 

Data from ATCA/ASKAP and EVLA could also be com- 
bined, but this makes little sense because the frequency cover- 
ages are similar. Data from ATCA/ASKAP and LOFAR are hard 
to combine because they cover diff'erent parts of the sky. 

The Global Magneto-Ionic Medium Surve y (GMIMS) 
planned with several single-dish telescopes (Wolle ben et al.l 
2010b) will provide a similarly wide frequency coverage as 
WSRT or GMRT, but with lower angular resolution. 

3. Wavelet-based RM synthesis 

The wavelet transform of the Faraday spectrum F((p) is given by 



00 



Wf(a, 0) = 




—00 



where i//((p) is the analyzing wavelet, a defines the scale, and 
(p is the Faraday depth of the wavelet center. We use as the 
analyzing wavelet the real- value "Mexican hat" i/^((p') = (1 - 
0'^)exp(-0'^/2). The coefficient wf gives the contribution of 
the corresponding structure of scale a located at position to 
the function F. Wavelet-based RM synthesis was introduced 
by Frick et al. ( 2010) and directly calculates the coefficients wf 
from the complex- valued pola rized intens i ties P (A^). The algo- 
rithm is described in detail by iFrick et aP (1201 ll) and allows us 
to combine the RM synthesis procedure and the wavelet filter- 
ing. With the additional assumption of symmetry of the radio 
sources, one can recover more accurately the Faraday spectrum 
F((p), particularly its complex part yielding the intrinsic polar- 
ization angles. 

We take advantage of the wavelet coefficient distribution as 
a multi- scale representation of the signal. The scale analysis 
is worth performing when the range of recognizable scales in 
Faraday space is sufficiently wide, in other words, the ratio of 
maximum to minimum wavelengths of the observations is con- 
siderably larger than unity 

mm mm 

Tabled shows that the EVLA, the ATCA, the SKA, and a com- 
bination of LOFAR, WSRT/GMRT, and EVLA fulfil this condi- 
tion. The SKA will provide the largest value of (/Imax/^min)^ in 
its phase 1, which will be further increased in phase 2. 
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Fig. 4. Wavelet planes wpia, 0) for the magnetic field with a Gaussian profile and in energy equipartition with cosmic rays (model 
1). The calculations are done for the frequency bands available for several telescopes. Left column, panels from above: LOFAR, 
EVLA (or ATCA), and ASKAP (POSSUM + FLASH); right column, panels from above: WSRT (or GMRT), ASKAP (POSSUM), 
and SKA phase L 



We note that a large value of (/Imax/^min)^ is also needed to 
obtain a small error Ai//o of the intrinsic polarization angle 

A^o = A0iL- V3/(|(|2i_i)), (5) 

min 

where S /N is the signal-to-noise ratio of a "source" at Faraday 
depth (p. 

4. Models 

We first consider a model region hosting a regular magnetic field 
with diff'erent profiles along the line of sight, as shown in Fig.[T] 



The maximum regular field strength of 2 fiG is typical of the 
Milky Way and nearby galaxies. The modeled region also con- 
tains thermal electrons (responsible for Faraday rotation) and 
cosmic-ray electrons (CRE, responsible for synchrotron emis- 
sion). Thus, the investigated region emits polarized synchrotron 
emission and rotates its polarization angles according to the 
Faraday eff'ect. 

The distributions of relativistic electrons Uc and thermal elec- 
trons n± are assumed to have Gaussian profiles in the radial 
direction within the plane, or with height perpendicular to the 
plane 

ndx) = CQxp[-(x- xof/hll (6) 
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Fig. 5. Wavelet planes wp{a, (p) for the magnetic field with Gaussian profile and the same scale-heights of thermal and CRE electrons 
(model 2). Panels are the same as in Fig.lH 



^thW = 0.03 cm ^ exp[-(x - x^f Ih^^. 

Here C is the maximum CRE number density, which is normal- 
ized to 1 in our models. Our models are applicable to galaxies 
observed at any inclination angle. If the scale-radius is, for exam- 
ple, ten times as large as the scale-height, the scale-height dom- 
inates the distribution for inclinations smaller than 84°, which is 
otherwise dominated by the scale-radius. 

We assume that the Gaussian scale-height (respectively the 
scale-radius) of the thermal electrons /zth is half that of the scale- 
height (respectively scale-radius) of the regular magnetic field 
/zb . The synchrotron scale-height /zgyn is derived according to 

(l//Zsyn)' = (l//^c)' + (2//ZB)'. 



For the CRE scale-height (resp. scale-radius) h^, four models are 
considered (Table 2): 

(1) /zc = ^hi\^ = h^l V2. This is the standard case expected 
for equipartition between the energy densities of cosmic rays and 
the regular magnetic field. The region of CRE is more extended 
than that of thermal electrons. The resulting /zsyn is the same as 
that of the thermal electrons. 

(2) he = h± = /^b/2. The regions of CRE and thermal elec- 
trons have identical extents and h^yn = V^T^ /zth- 

(3) he = /zth/ V2 = /zb/(2 V2). The region of CRE is less ex- 
tended than that of thermal electrons and h^yn = ^/2/5h±. This 
case is expected to occur at high frequencies and/or in regions 
with strong total magnetic fields where the energy losses of CRE 
are high and hence their propagation lengths are small. 
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(4) he = 2 hth = hB . The region of CRE is much more ex- 
tended than that of thermal electrons and h^yn = ^J4/3h±• This 
case is expected at low frequencies and/or in regions with weak 
total magnetic fields where CRE lifetimes and hence their prop- 
agation lengths are large. 

Changing the CRE scale-height leads to strong diff'erences in 
the shape of the Faraday spectrum. The first model results in a 
broad Faraday spectrum F with two strong "horns" at each end 
(Fig.O. These horns are the result of the two regions with very 
low density of thermal electrons on the far and the near side of 
the pathlength. One horn is correspondingly located at Faraday 
depth = (Fig. [T]). These horns disappear for the second 
model, in which both distributions have the same scale-height. 
However, the Faraday spectrum F still has sharp edges. The 



third model with a narrow CRE distribution has neither horns 
nor sharp edges. On the other hand, the extended CRE distribu- 
tion in model 4 leads to very strong horns. 

We conclude that the contrast between the amplitudes of the 
horns and of the broad part of the Faraday spectrum is controlled 
by the relative distributions of cosmic-ray electrons and thermal 
electrons. The larger the ratio hc/hth, the stronger the horns. No 
horns are visible for he < h±- 

In addition, we consider two cases with magnetic field rever- 
sals. The distributions of relativistic and thermal electrons are 
the same as in model 1, but the magnetic field includes one or 
two reversals (see Fig.[T] upper panel). Reversals along the line 
of sight can be the result of an axi- or bisymmetric spiral mag- 
netic field structure in the disk when a galaxy is observed almost 
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(rad m 



edge-on, or of an antisymmetric field structure with a reversal of 
the toroidal field component above and below the galactic plane 
when a galaxy is observed almost face-on. In the Milky Way, one 
field reversal is observed, which is locate d inside the solar radi us 
and extends over several kpc in azimuth (IVan Eck et al ]|20ni). 

The typical signature of the model with one reversal is the 
strong asymmetry of the Faraday spectrum F, called a "Faraday 
caustic" by Bell et al. (201 1), which becomes symmetric for two 
reversals, with two narrow features at the edges that are much 
stronger than in the basic case without reversals. 

Finally, having in mind objects such as the intergalac- 
tic medium in a galaxy cluster, we consider the case of a 
Kolmogorov-type turbulent field, without a regular field (Fig.©. 
The standard deviation ^RB^) = 5 jiG is typical of the intraclus- 
ter medium near the center of a radio-bright galaxy cluster. The 



field strength and both thermal and cosmic-ray electron densi- 
ties are assumed to have Gaussian distributions along the line of 
sight with scale-radii he = /zth = / V2. 

5. Results 

5.1. Recognizing regular fields and reversals in magnetized 
regions 

As noted above, the wavelet decomposition gives a graphical 
representation of a structure in the scale versus Faraday depth 
plane. We present in Figures l4lfT2l the modulus of the wavelet 
decomposition of the Faraday spectrum (Eq. lU, reconstructed 
from the range of provided by various radio telescopes, for 
the model examples described in the previous section. The hor- 
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Fig. 8. Wavelet plane for the turbulent field. Panels are as in Fig. |4] 



izontal axis in all panels gives the Faraday depth and the ver- 
tical axis shows the scale of the reconstructed structure (also in 
Faraday space). The color in the wavelet plane corresponds to 
the amplitude 

For simplicity, we assume for all wavelet transforms a con- 
stant polarized intensity of the emitting source and a constant 
wavelength response of the telescopes, to demonstrate the eff'ects 
of wavelength coverage. If the emitting source has a power-law 
frequency spectrum with a negative spectral index, the lower fre- 
quencies get a higher weight, so that the smaller scales in the 
wavelet plane are enhanced with respect to the larger scales. 
On the other hand, Faraday depolarization increases towards 
low frequencies, so that the frequency spectr um of polarized in- 
tensit y can reveal a positive spectral index (lArshakian & Beckl 
l20TTh . 

The wavelet plane for the first model (Gaussian profile of 
the regular magnetic field, he = V2/ze) is shown in Fig. IH One 
horn is located at = 0, while the location of the second horn 
gives the total Faraday depth of the region. The amplitude of the 
two horns on small scales is larger than that of the broad struc- 
ture, which we call "disk" in the following. We note that LOFAR 
can only recognize the two horns that cannot be distinguished 
from two point sources in Faraday space. Observing with the 
WSRT range picks up some of the medium- scale structures, but 
the small wavelength coverage and the large gap between the 
two frequency bands leads to high sidelobes that need to be re- 
moved with "RM clean". With the range of short wavelengths 
of ATCA, EVLA, and ASKAP (POSSUM) only the large scales 
are observable. Inclusion of FLASH indicates the presence of the 
two horns. However, only two broad structures are detected with 



ASKAP, which cannot be distinguished from two emitting and 
rotating regions with similar properties located along the line of 
sight. Only the SKA phase 1 can recognize all scales and is close 
to a perfect "Faraday telescope". Extending the frequency cov- 
erage of the SKA to 10 GHz (0.03 m) in phase 2 would help us 
to observe galaxies with large |A0|. 

The wavelet plane for the Gaussian profile of the regular 
magnetic field and = (model 2) is shown in Fig. [S] The 
sharp edges of F (see Fig.O generate horns in the wavelet plane 
as in the case of model 1 , but the amplitude of the horns is lower 
than that of the "disk". This is clearly seen with the wide wave- 
length coverage of the SKA phase 1. Observations with ASKAP 
are already helpful in this case because the scale of the transi- 
tion between disk and horns lies in the range of scales covered 
by POSSUM + FLASH (but not with POSSUM data alone). The 
other telescopes are unable to distinguish between models 1 and 
2. 

The wavelet plane for a distribution of a regular field with 
one reversal is shown in Fig. |6l The dominating horn is eas- 
ily detectable with LOFAR, WSRT/GMRT, and the SKA phase 
1 . However, observations at long wavelengths alone cannot dis- 
tinguish this structure from a single point source in Faraday 
space. Other than for a point source, the structure of the horn 
is asymmetric, which can be recognized by the inclusion of data 
at higher frequencies, such as with the SKA or by combining 
POSSUM -h FLASH. On the other hand, the response from high 
frequencies alone (ATCA, EVLA, and POSSUM) is a somewhat 
asymmetric disk that cannot be distinguished from e.g. asym- 
metric Gaussian profiles. The second horn is weaker by a factor 
of ten and is hardly detectable. Observations at several telescope 
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Fig. 9. Wavelet plane for the turbulent field superimposed on the regular field. Panels are as in Fig. |4l 



pointings or complete mapping of the source is needed to deter- 
mine the structure of sources and the extent of field reversals in 
the resulting data cube. 

The wavelet plane for a distribution of a regular field with 
two reversals, which may exist in spiral galaxies, is shown in 
Fig. 13 The horn to the left is shifted to negative cp because of the 
negative components of the magnetic field along the line of sight 
(Fig.[T]). The "disk" at large scales is very weak compared to the 
horns, even weaker than in Fig. |6l and is hard to detect at short 
wavelengths. At long wavelengths (LOFAR and WSRT/GMRT), 
the two horns cannot be distinguished from two point sources 
(Faraday screens). The asymmetric nature of the two horns is 
visible at intermediate scales and can be recognized from the 
shift of the maximum at diff'erent scales. 



Double features in the Faraday spectrum were found in the 
central regions of several spi ral galaxies from the WSRT SINGS 
survey at 1300-1763 MHz ( Heald et al.ll2009l) . As the WSRT is 
sensitive only to large-scale structures in Faraday space, these 
features cannot be interpreted as field reversals and are probably 
distinct emitting and rotating regions along the line of sight, e.g. 
a regular fie ld component assoc iated with the nuclear region as 
proposed by lHeald et al.l (l2009l) . Observations at lower frequen- 
cies are needed to achieve higher resolution in Faraday space. 

5.2. Turbulent fields 

The signatures of turbulent fields are many components on small 
scales, which we call the "Faraday forest", can only be observed 
at long wavelengths (Fig. [8]). In addition, owing to the tail of the 
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Fig. 10. Wavelet plane for the turbulent field superimposed on the regular field with one reversal. Panels are as in Fig. |4] 



continuous turbulence spectrum on large scales, a weaker struc- 
ture is visible on yet larger scales ofa^S rad m"^. This structure 
is less extended in (p than the "disk" in the model with turbulent 
+ regular fields (Fig. [9]). However, observations at only high fre- 
quencies cannot distinguish the extended structures generated by 
either the tail of the turbulence spectrum or a regular field. 

In real observations, the emission is smoothed over the tele- 
scope beam and distorted by instrumental noise. The "Faraday 
forest" is most clearly visible when the spatial scale correspond- 
ing to the beam is similar to the turbulence scale or smaller. For 
a larger beam, the "Faraday forest" is less visible and causes de- 
polarization. 

Instrumental noise is spread over the whole Faraday spec- 
trum and hence is easy to distinguish from the "Faraday forest". 



even if the amplitude of the "forest" is similar to that of noise 
( Frick et al . 2011). 

The case of a Gaussian regular field superimposed on a 
Kolmogorov-type turbulent field is shown in Fig.|9l The standard 
deviation ^/{Bj} is assumed to be three times larger than the reg- 
ular field strength, which is approximately valid for spiral galax- 
ies. In contrast to Fig. [8l the Faraday forest extends only over a 
limited range in 0. On intermediate scales (WSRT/GMRT), the 
disk becomes asymmetric, but no clear signatures of the turbu- 
lent field are visible. 

The corresponding result for a regular field with one reversal 
-I- turbulent field is shown in Fig. [TOl The components of the 
"Faraday forest" are stronger at one edge of their distribution, as 
can easily be recognized at long wavelengths. 
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Fig. 11. RM synthesis applied to data from a combination of data from LOFAR and EVLA. Panels in the left column from above: 
box-like magnetic field, one reversal, and box-like field + turbulent field; right column from above: Gaussian-like distribution, two 
reversals, and regular field with one reversal + turbulent field. 



5.3. Recognizing magnetic structures by combining data 
from several telescopes 

Finally, we combine the wavelength ranges of LOFAR and 
EVLA and apply RM synthesis for our diff'erent field models 
(Fig.[TT]). The lack of intermediate wavelengths leads to a sepa- 
ration between small and large scales in Faraday space. Box- like, 
Gaussian, and double-reversal models can be distinguished from 
the relative amplitudes of the horns and the disk. 

By combining data from LOFAR, WSRT (or GMRT), and 
EVLA (Fig. [121), structures at intermediate scales can also be 
recognized and yield an almost complete picture of the magnetic 
field distribution, comparable to the SKA. 

6. Conclusion and discussion 

Our conclusions are as follows: 



1 . We have found that a reliable recognition of magnetic field 
structures in either the interstellar medium of spiral galaxies or 
the intracluster medium of galaxy clusters requires us to apply 
RM synthesis to spectropolarimetric data cubes observed with 
high angular resolution and over a wide frequency range, from 
about 100 MHz to several GHz. Such a wide frequency cover- 
age provides high resolution in Faraday space (as determined by 

as well as recognition of a wide range of scales in Faraday 
space (as determined by (/Imax/^min)^)- 

2. The combination of data from the POSSUM and FLASH 
surveys (both planned with the ASKAP telescope) would im- 
prove the recognition of structures on intermediate scales. 

3. The combination of data from the present telescopes 
LOFAR and EVLA appears to provide a promising means of 
recognizing magnetic structures on all scales. 

4. The combination of WSRT (or GMRT) data with those 
from LOFAR and EVLA would fill the gap between the LOFAR 



R. Beck et aL: Recognizing magnetic structures with RM synthesis 13 



0.00 0.02 0.04 0.06 0.08 0.00 0.01 0.02 0.03 0.04 0.05 0.06 




(p (rad m ^) ^ (rad m ^) 

Fig. 12. RM synthesis for a signal from a combination of data from LOFAR, WSRT (or GMRT), and EVLA. Panels are as in Fig.fTTl 



and EVLA frequency ranges and hence can also recognize inter- 
mediate scales in Faraday space, which is helpful for measuring 
magnetic structures. 

5. The detection of two "horns" on small and intermediate 
scales plus an extended "disk" on large scales in the wavelet 
plane of galaxies indicates that the scale-height (or scale-radius) 
of cosmic-ray electrons (CRE) is larger than that of thermal elec- 
trons, as expected especially at low frequencies. The amplitude 
of these "horns" relative to that of the "disk" allows the deter- 
mination of the ratio of the scale-heights (or the scale-radii) of 
both electron populations. To distinguish these horns from point 
sources in Faraday space, detection of intermediate scales would 
be required. 

6. The detection of an extended "disk" without "horns" in 
the wavelet plane of galaxies indicates that the scale-height (or 
scale-radius) of CRE is smaller than that of thermal electrons, 
owing to e.g. the strong energy losses of CRE, which are ex- 
pected at high frequencies. 



7. Recognition of field reversals in spiral galaxies needs de- 
tection of structures on small and large scales, e.g. by combining 
data from LOFAR -h WSRT/GMRT EVLA, or with the SKA, 
which is close to a perfect "Faraday telescope". 

8. Turbulent fields in galaxies or the intracluster medium can 
be recognized on small scales as a "Faraday forest" of many 
components, observable with high angular resolution and at long 
wavelengths, e.g. with LOFAR or the SKA. 

9. The single-dish, all- sky polarization survey GMIMS will 
provide excellent resolution and scale recognition in Faraday 
space. The low angular resolution allows us to investigate the 
structure of the magnetized medium in the Milky Way. 

10. For simplicity, the eff'ects of instrumental noise and vari- 
ations in both angular resolution and polarized intensity with fre- 
quency have been neglected in this paper. In practice, the combi- 
nation of data from diff'erent telescopes suff'ers from varying an- 
gular resolutions and, in the case of synthesis telescopes, from 
diff'erent distributions of baselines (uv coverage). Furthermore, 
both a non-zero spectral index of polarized intensity and the sig- 
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nal averaging within the telescope beam modify the visibility of 
scales in Faraday space. Understanding these effects will require 
more detailed modeling. 

1 1 . Mapping of the source is required to help us recognize 
the magnetic structures in the 3-D data cubes (spatial coordi- 
nates -\- Faraday depth coordinate). As the next step, the wavelet 
analysis will be extended to analyzing Faraday data cubes. 

12. We propose to call the data cubes generated by RM 
synthesis "PPF (position-position-Faraday depth) cubes", to 
be analogous to the PPV (position-position-velocity) cubes in 
spectroscopy. 
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Appendix A: 

Fig. lA. II shows the wavelet decompositions of the various mod- 
els for complete frequency sampling, as will nearly be reached 
by the SKA phase 2. The result is almost the same as SKA phase 
1, except for very large scales a > lOOradm"^. 



R. Beck et aL: Recognizing magnetic structures with RM synthesis 15 



0.00 0.02 0.04 0.06 0.08 0.10 0.00 0.01 0.02 0.03 0.04 0.05 0.06 




(rad m ^) ^ (rad m ^) 



Fig. A.l. Wavelet plane of the models considered in Sec. 5 but for complete coverage in frequency. Top: Regular field of model 1 
(compare with Fig.|4]) and of model 2 (compare with Fig.O. Middle: Regular field with one reversal (compare with Fig.|6} and with 
two reversals (compare with Fig. |7]l. Bottom: Turbulent field superimposed on the regular field of model 1 (compare with Fig.O 
and turbulent field superimposed on the regular field with one reversal (compare with Fig. [TOl l. 



